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Abstract 

Classical Lie group theory provides a universal tool for calculating symmetry 
groups for systems of differential equations. However Lie's method is not as much 
effective in the case of integral or integro-differential equations as well as in the case 
of infinite systems of differential equations. 

This paper is aimed to survey the modern approaches to symmetries of integro- 
differential equations. As an illustration, an infinite symmetry Lie algebra is calculated 
for a system of integro-differential equations, namely the well-known Benney equa- 
tions. The crucial idea is to look for symmetry generators in the form of canonical 
Lie-Backlund operators. 

1 Introduction 

The major obstacle for the application of Lie's infinitesimal techniques to integro-differential 
equations or infinite systems of differential equations is that the frames (see, e.g. or 
0) of these equations are not locally defined in the space of differential functions. In 
consequence, the crucial idea of splitting of determining equations into over-determined 
systems, commonly used in the classical Lie group analysis, fails. 

1.1 Different forms of the Benney equations 

The Benney equations referred to by the name of the author of a pioneering work 
appear in long wavelength hydrodynamics of an ideal incompressible fluid of a finite depth 
in a gravitational field. From the group theoretical point of view they are of particular 
interest due to the existence of an infinite set of conservation laws obtained in Q. The 



latter property of the Benney equations emphasizes their significance that goes far beyond 
an interesting example of an integrable system of hydrodynamic equations. 

In practice, the Benney equations are used in various representation. One of them is 
the kinetic Benney equation (a kinetic equation with a self-consistent field): 

+00 

ft + vf, - Alf, = 0, A\t, x)= j fit, X, v)dv. (1) 

—00 

This equation appears as a unique representative of a set of hierarchy of kinetic equations of 
Vlasov-type A detailed study of its group properties will lead to better understanding 
of the symmetry properties of kinetic equations of collisionless plasma, viz. the Vlasov- 
Maxwell equations (see also Appendix 1) that have both theoretical and practical interest, 
e.g. while dealing with controlled nuclear fusion programme. 

Another form of the Benney equations is an infinite set of coupled equations 

4 + + iAlA^-^ = 0, f > (2) 

for a countable set of functions A* of two independent variables, time t and the spatial 
coordinate x. In terms of hydrodynamics these functions appear as averaged values (with 
respect to the depth) of integer powers i > of the horizontal component of the liquid 
flow velocity. The corresponding integrals that describe this averaging are taken over the 
vertical coordinate in the limits from the flat bottom up to the free liquid surface. Solu- 
tions, Hamiltonian structure and conservation laws for the equations (|^ were discussed 
in details in |^, |6| . 

^From the kinetic point of view the system (^) can be treated as a system of equations 
for moments of the distribution function / that obeys the kinetic Benney equation (||) 

+00 

A\t,x) = j v'fdv, i>0. (3) 

— CO 

This fact with the explicit formulation of the Benney equation (|l|) was first stated inde- 
pendently in 0, g. The La grangian change of the Euler velocity v, 

v = V{t,x,u) (4) 

yields one more representation for Benney equations (||): 

ft + Vf, = , Vt + VV^ = -Al , AO(t, x) = J Vufit, X, u)du. (5) 
The equations (|5|) are readily converted into the hydrodynamic-type form 

nt + {nV)^ = Q, Vt + VV^ = -Al, A° = J n{t,x,u)du, (6) 
if one employs the "density" n depending on the Lagrangian velocity u: 

n = f{t,x,u)Vu. (7) 
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Using the form (|6|) of the Benney equations an infinite set of conservation laws were 
constructed in with the densities regarded as functions of the Lagrangian velocity u. 

Furthermore, we will rewrite the integro-differential Benney equations ([l|) in the form 
of differential equations by introducing the following nonlocal variables: 

V +00 

g= j fdv, h= J fdv. (8) 

—00 V 

In terms of the latter variables the equations (|^) are written as 

ft + vf^-{ga: + K)fv = 0, gv = f , K = -f- (9) 

The knowledge of the complete Lie-Backlund symmetry for the Benney equations in 
different representations ([l|), (||), @, (P) and (P) can clarify the question of structure of 
solutions and conservation laws for these equations. This statement is partially confirmed 
by the fact that one of the main results of the works |^ ^, namely the higher order 
Benney equations, can be re-formulated in terms of the first order Lie-Backlund group, 
admitted by the system (||) . Unfortunately, the complete description of the Lie-Backlund 
symmetry for the equations (||) is not available in the literature. The goal of this paper is 
to contribute to this problem by calculating an infinite (countable) part of the Lie point 
symmetries of the moment equations @. 

2 Generalities 

Here, we briefly discuss different known approaches to calculating symmetry groups for 
integro-differential equations. Loosely speaking, these approaches can be divided into two 
large groups: indirect and direct methods. 

Algorithms of the first group rest on the possibility to replace in any way input nonlocal 
(integro-differential) equations by a system of differential equations. Then the resulting 
system of differential equations is analyzed using standard methods of a classical Lie group 
analysis. Here we point on two different ways of reducing nonlocal equations to differential 
ones. 

2.1 Indirect methods 

A. Method of moments 

In this approach, the system of basic integro-differential equations for a function / (e.g., 
the kinetic equation (0)), that usually contains nonlocal terms depending on moments (|3|) 
of this function, is reduced to an infinite system of differential equations for these moments 
(in our case this is the system (^)). The admitted symmetry group is then calculated us- 
ing the traditional methods of Lie group analysis for any finite subsystem of k equations 
of this general infinite system. Then an intersection of all admitted groups is defined 
and a transition to the limiting case A: — > 00 is fulfilled. The resultant algebra of group 
generators thus obtained is used to reconstruct the algebra with the original function / 
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directly involved. This procedure usually employs the explicit form of finite group trans- 
formations for moments (^) and the relations between A* and the function /. The last 
step is not trivial in any case as there may exist different representation for the equations 
for the moments and the resultant group depends on the form of this representation. 
Hence, the transition from the algebra of group generators for the moments representation 
back to the symmetry of original equations formulated in terms of input functions form 
a special problem which we are going to discuss in the section The above described 
algorithm was realized to calculate Lie point symmetry group for Vlasov-Maxwell equa- 



tions in plasma theory ||9[-|11] and for Benney, Vlasov-type and Boltzmann-type kinetic 
equations p|, |]. 



B. Method of boundary- differential equations 

The following method was developed in ^] on basis of the concept of covering and ap- 
plied to a coagulation kinetic equation. In this case each definite integral in the equation is 
replaced by the corresponding difference of values of the antiderivative on a boundary sets. 
The integro-differential equation takes the boundary (or functional) differential form. As 



shown in [14, 15 1 a geometric theory of boundary differential equations can be constructed 
in just the same way as the analogous theory of differential equations thus allowing to 
define and compute not only classical but higher symmetries as well. It is noted that 
the elimination of integrals by introducing potentials depends on the choice of potential 
variables and can be executed in different ways. Hence the result of group calculation and 
its dimension is influenced by the form of potential variables involved. 

2.2 Direct methods 



Direct methods of finding symmetries were developed in [ 16 1- 1211] (^^^ ^^^^ KSl) ^^'^ ^1 
and applied to find symmetries of kinetic Boltzmann equation, the equations of motion 
of viscoelastic medium and Vlasov-Maxwell equations of plasma theory. To extend the 
classical Lie algorithm to integro-differential equations it appears necessary to resolve 
several problems. First, one should define the local one-parameter transformation group 
G for the nonlocal (integro-differential) equations and formulate the invariance criteria 
that lead to determining equations. 

Let an integro-differential equation under consideration be expressed as a zero equality 
for some functional (here we indicate only one argument for a function /) 

F[f{x)]=0, (10) 

and let G be a local one-parameter group that transforms / to f{x), 

f (x) = f -\- + 0(0"^) , x = x, (11) 

Here we use the canonical group representation hence independent variables x do not 



vary. Then the local group G of point transformations (11) is called a symmetry group of 



integro-differential equations ( p!o| ) iff for any a the function F does not vary [17|, 

f[/(x)1=0. (12) 
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Differentiating (12) with respect to group parameter a and assuming a — > gives the deter- 
mining equations. In contrast to the case of input differential equations these determining 
equations are in general also integro-differential. 

The invariance criterion for F with respect to the admitted group can be expressed in 
an infinitesimal form using the canonical group operator y, 

YF 



= 0, where Y = I dysi {y) -—— . (13) 
F=o J 6f (y) 

Here with the goal to generalize the action of a canonical group operator not only on 
differential functions but on functionals as well we use variational differentiation in the 



definition of y [23, 24 1. One can verify by direct calculation that the action of Y on 
any differential function and its derivatives, e.g., / and fx, ■ ■ ■ produces the usual result: 
Yf = ae, Yfx = Dx{se) and so on. Hence, if F describe usual differential equations 



then formulas (13) lead to standard local determining equations, while for F having the 
form of integro-differential equations formulas ( |l3|) can be treated as nonlocal determining 
equations as they depend both on local and nonlocal variables. 

In order to find solutions of determining equations one can use different approaches, 
e.g. expanding coordinates of group generator into formal power series and equating coeffi- 
cients of various powers [|l^, 0]. However there exists a more traditional way. As we treat 
local and nonlocal variables in determining equations as independent it is possible to sep- 
arate these equations into local and nonlocal. The procedure of solving local determining 
equations is fulfilled in a standard way using Lie algorithm based on splitting the system 
of over-determined equations with respect to local variables and their derivatives. As a 
result we get expressions for coordinates of group generator that define the so-called group 
of intermediate symmetry [^3|. In the similar manner the solution of nonlocal determin- 
ing equations is fulfilled using the information borrowed from an intermediate symmetry 
and by "variational" splitting of nonlocal determining equations using the procedure of 
variational differentiation. Therefore, the algorithm of finding symmetries of nonlocal 
equations appears as an algorithmic procedure that consists of a sequence of several steps: 
a) defining the set of local group variables, b) constructing determining equations on basis 
of the infinitesimal criterion of invariance, that employs the generalization of the definition 
of the canonical operator, c) separating determining equations into local and nonlocal, d) 
solving local determining equations using a standard Lie algorithm, e) solving nonlocal 
determining equations using the procedure of variational differentiation. 

In the next sections, the above methods are applied to the Benney equations. 

3 Lie subgroup and Lie-Backlund group: statement of the 
problem 

3.1 Lie subgroup: direct method of calculation 

A Lie subgroup, admitted by the kinetic Benney equation (|^) in the space of four group 
variables 

t, X, V, f (14) 
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is defined by five basic infinitesimal operators 

= dt, X2 = dx, X3 = tdx + dv, 

(15) 

X4 = tdt - vdy - fdf, X5 = xdx + vdy + fdf . 

With the less computation difficulties this group can be obtained using the approach 
developed in |2^, Q in application to group analysis of Vlasov-Maxwell equations in 
plasma theory (see also Chapter 16 in Handbook ||2^ ). 

Here we demonstrate the application of the general scheme to Benney equations ([l|). 
The merit of this direct approach is the possibility to present the criterion of invariance 
with respect to local one-parameter point group transformations in a standard infinitesimal 
form. For the Benney equations the point symmetry group generator has the form 

X = ^dt + i^dx + fd, + Tj^df + rj^d^o , (16) 

where the coordinates ^ and ry depend on t, x, v, f and A^. 
In the canonical form this operator is written: 

Y = si^df + 86^5^0 , (17) 

where 

^' = ri'- eft - eu - , ^' = v'- - eAi , 

and its action on any function or functional should be understood in generalized sense as 
in (13). Applying the canonical group operator to the joint system of basic equations @i 
and one more "evident" equation, that expresses the fact that the moment of / does not 
depend upon v 

^2 = 0, (18) 
gives the system of determining equations 

A(ae^) -^S^^(se^) -^x(«^)^ = 0, (se^) = ; (19) 



(B^dv, (20) 



which should be solved in view of the the complete set of basic equations (jl]), (|T8[). 
In solving determining equations ( |T9|) group variables 

{f,A^,fx,Al,f,,...} (21) 

are treated as independent ones. This assumption separates determining equations into 
local equations (|^) and nonlocal equation (|20|). Local determining equations are solved 
in a standard way using the computational algorithm of Lie group analysis. Then the 
functions ^ and r/ thus obtained define the so-called intermediate symmetry (Esf 



(22) 



= riHf) , = lit) + {2a - e)A' - xPu - ^ Clu ■ 
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Here C^{t), [5{t), 7(t) and rj^{f) are arbitrary functions of their arguments and a and u 
are constants. 

Now let us turn to the solution of the nonlocal determining equation (^0|) that we 
rewrite in the following form 



r] 



eA', - =J{v'- eft - - eu) dv . 



(23) 



As in the case of the local determining equations (0), the latter should be solved in view 
of the original equations (||), (|l^). Hence, calculating the derivatives A^, and ft from 
the basic equation (|l|) and inserting them into (|2^) and in view of the above expressions 
for coordinates ^ and rj we obtain the following nonlocal determining equation 



J L\f)-(a-^(Af 



dv = -f- x(3tt - — iltt 



(24) 



As any determining equation, ( p4D is the equality with respect to all group variables that 
appear in this equation. Therefore differentiating it with respect to any group variables 
also leads to equalities. Hence, nonlocal determining equation can be split with respect to 
independent group variable / using the variational differentiation. Since the right-hand 
part of (24) does not depend on /, it reduces after differentiation and the remaining terms 
are written as 



5f{v') 



n\f) 



a 



\e]f 



dv = 0. 



(25) 



It is essential that the nonlocal determining equation ( |23| ) should be solved simultaneously 
with its differential consequence, i.e. any solution of (|25| ) must appear as the solution of 
(pi). Introducing the variational derivative {5/5 f{y')) inside the integral over v 



+00 



5f{v 



Sf{v^ 



and eliminating integration over v with the help of the Dirac delta-function, 

Sf{v) 



5f{v') 



5{v - v') , 



one comes to the first order differential equation for rj^ 



0, 



which gives the linear dependence of rj^ upon /, 

1 



(26) 
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with some constant C. Substituting this result back into (p^) yields the zero-value of this 
constant provided integral has finite value. As r/^ does not depend on t then differentiating 
(p6| ) with respect to t gives .^^^ = 0. Differentiating (^) with respect to x gives two more 
equations, 

7 = Pu = 0. 



Solving these equations gives final expressions for coordinates and ij: 

^^ = c^+cH, = c^ + ch + e'x, f = + {-c-^ + c^)v , 



(27) 



These coordinates give rise to the five-dimensional Lie algebra with generators given by 
(p^). One can see that terms proportional to which come from se^ are omitted in 
X4, X5 as they appear as the result of prolongation [^] of these operators on a nonlocal 
variable A^. This procedure for any of operators X4, X5 is described in a concise form in 
the next section. 

3.2 Extension of Lie point symmetry generators to nonlocal variables 

To fulfill the procedure of prolongation of the Lie point symmetry generators one should 
first rewrite the operator, say X^, in a canonical form 

Y5 = a,^df , = {-xU - vh + f) ■ (28) 

Then formally prolong this operator on the nonlocal variable A^ 

¥5 + 3}^ = si^df + si^dAo . (29) 

The integral relation between ae^ and ae^ is obtained by applying the generator (^9|) to 
the second equation in that is treated here as the definition of A^. This relation here 
coincides with nonlocal determining equation (pc|). Substituting ae^ from (^) in (20) and 



calculating integrals obtained (integrating by paths) gives the desired coordinate se^ 

= {2A^ - xAl) . (30) 

Inserting this into (^) and returning back to the non-canonical representation we get the 
following generator 

X5 = xd^ + vd^ + fdf + 2A^dAo , (31) 

that correlates with the result (p7|). 

Prolongation of infinitesimal operators (ITsI ) on nonlocal variables (^) extends the set 
of group variables (|l4|) up to a countable set 

t, x, V, f, A^, A\ ... . (32) 

In the latter case infinitesimal operators ( [l5| ) rewritten in the canonical form |]l|, Sec. 8.4.2.] 
and restricted on the sub-manifold 

^, X^ A ^ • • • ; 

A\ .... (33) 
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are given by the following expressions 



i=0 1=0 i=0 

oo oo 

i=0 j=0 



(34) 



It can be easily checked that infinitesimal operators ( |34D are admitted by Benney equations 
(|2|) and it goes without saying that they directly result from the group analysis of Benney 
equations (|2|). Just in this way (i.e., using the method of moments) infinitesimal operators 
(|l5| ) were first obtained in by using non-canonical form of infinitesimal operators (| 
with the subsequent passage to the representation (p^) in the space of variables (p2D. 

3.3 Incompleteness of the point group: statement of the problem 



It is evident, however, that the subgroup (l34|) does not exhaust the complete group sym- 
metry of Benney equations (^). The incompleteness of the result ( |3^ ) is obvious form 
many points of view. Here we shall only point on the nonconformity of finite dimension 
of the algebra ( |3^ to the infinite set of conservation laws for Benney equations, and on 
the infinite extension of the point symmetry group for Benney equations in the form of 



(|5|), (y) with Lagrangian velocity (see also in the context Chapter 16 in |25| where this 



extension was outlined for Vlasov kinetic equation in plasma theory). Here of principle 



significance for us is the following statement: the group (34) is incomplete not only from 
the standpoint of Lie-Bdcklund symmetry for Benney equations but also from the stand- 
point of the Lie point symmetry. The validity of the statement can be proved by direct 
solving of determining equations for the first order Lie-Backlund group (contact group, 
that is not reduced to point one) 

Dt{si')+D^{cs'+^) + iA'-^D^{cS^) + iAlcs'-^ = 0, i>0, (35) 

where coordinates of canonical operator 

00 

X = J2 ^'d^i , (36) 

i=0 

depend upon the countered set of group variables 

t, x; A^ , . . . , A-^ , . . . ; A^, . . . , A^, . . . ; j > . (37) 

To prove the above statement one can consider only partial solutions of determining equa- 
tions (H) 

si' = r,\A^,...,A^,...); i,j>0, (38) 

that depend upon moments A^ , j > 0, and does not depend upon t, x. It appears that 
thanks to these infinitesimal operators (|36|), (^) an infinite extension of the group (| 
takes place. Now the problem is to find these operators. 
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4 Determining equations and their solution 

4.1 General form of the determining equations 

Before proceeding further we write determining equations of first-order Lie-Backlund 
group, admitted by a more general (as compared to (|2|)) infinite system of coupling equa- 
tions for functions A'^{t,x) with the arbitrary element (p{A^) 



Ai + ^^+1 + iA'-^ 



ip{A'^) =0, i>0. (39) 



For the coordinates ae* of canonical infinitesimal operator (|^) the following chains of 
determining equations are valid which result from splitting (^9[) with respect to second 
derivatives: 

oo 
j=0 

seji + iy^iA^-ia^o , = ee\^, ; i > , j > 1 , 



aej + ae^+^ + iipiA^^'^sil + A° («(/9i8e*~^ + iip2A^''^si 

oo 

+ J2 [iv^iA'-'Ais,% - {Ai+' + MAlA^-') ^% + ^^ae 

oo 

jAl ((^i^r 1 + ^2AlA^-^) 4, = , i > 



(40) 



i+l 
Ai 



3=0 

Here (pi and 992 are the first and the second derivatives of the function ip with respect to 
its argument. From the various standpoints at list three distinct values of the function 
are specified. In case 'p^{AP) = AP we come to kinetic Benney equations (|l]), whereas 
for if = a^A^)"^ extension of the admitted point group takes place thanks to projective 
transformations in t,x-plane (see [^). For tp = ahiAP the corresponding kinetic equation 

A^ Y 
ft + vU-a^U = Q, A'=jdvf, (41) 

—00 

that gives rise to the discussed system of equations for moments, is of special interest in 
plasma theory. It appears as the equation for the distribution function of plasma ions, 
while electrons obey the Boltzmann distribution. More complicated dependencies of p>{A^) 
upon A'> can also be of interest in plasma physics for non-Boltzmann distribution functions 



for hot electrons. The equation (41) was studied in details in ||2] 

For the Benney equations (^) the determining equations (|40| ) are rewritten in the 
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(42) 



following form 

oo 

asjii - as^, + iA'-^sP^+, = , i > , j > . 

(oo oo 

oo oo oo 

+ z^i-i E + E - E = , i > . 

i=o j=o i=o 

4.2 Solution of the determining equations 

Under conditions (38) the determining equations (^) are splitted and reduced to two 
infinite chains of equalities, namely one-dimensional (vector) and two-dimensional (tensor): 

oo 

^aV - E + ^^^"Sao + = , ^ > ; 

i=o (43) 

rf+\, - ri% + iA'''v%+i = , i > , > . 

The apparent difficulty in analytical solving of the given system of determining equations 
( |43| ) is due to a "nonlocal" nature of the second term in the vector chain in the form of an 
infinite sum with respect to index j > 0. The measure of this non-locality is characterized 
by a number of nonzero components of tensor r/j. But in fact in case of an overdetermined 
system ([43|) we obtain a finite upper value of the summation index j < oo, which depends 
upon the other index i of this tensor. In order to be sure that it is true we shall first 
consider the system of determining equations ( |4^ ) in two particular cases, namely for 
i = and i = 1. In case z = we have two coupled determining equations just for two 
coordinates rj^ and rj^ of the desired infinitesimal operator (|36D, (| 



r/]io -Ei^'"'<^ =0' v\k+i=V%, k>0. (44) 

3=0 

The second determining equation in ( ^4| ) enables to eliminate the coordinate rj^ from the 
first determining equation and obtain as a consequence of the system (44) the following 
isolated scalar determining equation for the coordinate t]^ only 

oo 

^AO-Ei^'"'^A.+i=0. (45) 

j=0 

For i = 1 the system (p3|) yields coupled equations for three coordinates rj^, rj^ and 77^ 

00 

r/io=Ej^'"'^i^-(^° + ^^Ao), V%+^=V\^-A'^r^%+^, k>0. (46) 

j=0 
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Compatibility conditions for determining equations 

??^o^fe+i = f?^fc+i^o, A;>0, (47) 



enables to eliminate the coordinate rf from ( ]46| ) and obtain one more closed vector deter- 
mining equation for ry^ (more precisely, the determining equation that contains the first 
and the second derivatives of r]^ with respect to A^): 

oo 

^V-Ao - E - + 2)'7a/=+2 = , > . (48) 
i=o 

Differentiating the scalar equality (^) by moments yields one more vector corollary 
for the coordinate r/^: 

oo 

^AOAfe - E j^^"^v\,+iA^ ~{k+ l)ri\^+2 = , > . (49) 



From the compatibility conditions for two determining equations ( |4^ ) and ( p9D it follows 
that the coordinate r/^ does not depend upon the moments A^ which are higher than A^ 

r]\.+2 = , i > . (50) 



This formula arises as a result of mutual subtraction of determining equations ( |48| ) and 
(]49| ) in view of the equality of terms containing summation over index i > 0, that is the 
corollary of the second determining equation in ( ^4| ) 

v\j+iA^ = v\jAk+i ' i, A; > . (51) 

The symmetry of the second derivatives of r/^ given by ( |5ll ) results from the tensor form 
of a compatibility condition for two vector determining equations from (^), which differ 
only in "sounding" index 

r/^fc = , V% = , j, k>0. (52) 

The result ( [50| ) is a milestone on the way of solving the system of determining equations 
(|43|). Indeed, in view of (|50| ) the second determining equation of the system (^) yields 
the requirement of independence of the coordinate t]^ upon all higher moments, that differ 
from A^ 

V%+i=0, i>0, (53) 
and with ( |53| ) in mind the tensor chain from (43) is simplified in such a way 

=r/>, i>0, k>0, (54) 



that enables to modify the equality (|50| ) in the sense that any coordinate r/* of the desired 
operator (|3^, ( |38| ) does not depend upon any moments A^ , higher than A^ 

V%+,+i = , i>0, k>0. (55) 
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The last equality sets a finite upper limit i > j to the summation index j > in the second 
(nonlocal) term of the left hand side of the vector determining equation of the system ( p3|) 



^aV - E + iA'-'vAO + ir]'-' = , i > . (56) 

3=0 

The use of one of the equalities ( |5^ ) 

V%=ri%, i>0 (57) 
enables to rewrite the chain of determining equations ([5^) in even a more simple way 



j=0 



(58) 



whereas one of the corollaries of determining equations (54), that results for i = A; + 1 in 
combination with (53) 

r?^V=0, i>0, (59) 
lowers the upper value of the summation index j > in (^) by unit 



i-2 



Va' - E J^^'Wa. + iv'-' = , i>0. 

j=0 



(60) 



Collecting the arising determining equations (p^), ( [5^ ) and (|59| ), ( pOD we arrive to a much 
more simplified (but equivalent) formulation of the system (|43|), which really is integrated 
below 



i-2 



i=o 



r/^V+i =^V-, r/^.+fc =0, i>0, /c>0. 



(61) 



The last of the four equalities in ( |6l| ) is strengthened in comparison with (|5^) thanks to 
the condition 

V%=0, i>0, (62) 



which can be obtained as follows. For example, the first two equalities (62) for i = and 
i = 1 result as a corollary of (57) and compatibility conditions 



(63) 



for the first derivatives of the coordinate r]^ with respect to moments and (see ( |60[ ) 
for i = 1 and i = 2) 

r,% = -27?i , r,% = 7j\o = -rf . (64) 



After that, the validity of the remaining equalities (p3) becomes obvious thanks to (|5^). 
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Before proceeding to enumerating all solutions of the system of determining equations 
(|6l|), we present here yet another form of the chain ( |60|) 

i-2 

rj^,^-Y^jA^-^rj%^+i^^-^ = 0, i>0. (65) 

j=0 

This form can be employed to clarify the general structure of these solutions on basis of 
the corresponding generating functions. 



4.3 Discussion of the solution of the determining equations 

The integrability procedure in itself for determining equations (|6l|) is of no difficulties. 
For example the first six coordinates ry* (0 < i < 5) of the desired infinitesimal operator 
(p6|), ( p8[ ) are given by the following formulas for the general solutions of determining 
equations (|6l]) that depend upon six arbitrary constants (0 < j < 5) and are described 
by polynomials in moments 



1 aO 



-A^ + {A 



Q\2 



(66) 



- AC^A^ - SC^A^ + [-2^2 + 3(A°)2] +C° (-A^ + 2A°A^) . 



It appears that the polynomial dependence of any solution rf of determining equations 
(|6l| ) upon moments A^ is a general property of components of the vector rf for any i > 0. 
The example (|6^) demonstrates that the procedure of obtaining solutions of determining 
equations (|6l| ) is reduced to their enumeration. To be concrete, we assume the following 
scheme of indicating of the k-ih. basic solution of determining equations ( |6l[ ) for the 
coordinate rf: 



Vk 




k, i>k + 2; i, k > . 



(67) 



In the solutions (66) this scheme demands quit definite choice of values of integration 
constants in the form of Kronecker symbols 



CJ 



"jjk ; 



j, > . 



(68) 



The last of the four equalities for rf^ in (|67] ) (in square brackets) indicates the homogeneity 
degree (i — k) of the polynomial "tail" of the solution rf for z > A; + 2 in accordance with 
the attributed to any of the moments A^ of the order i the homogeneity degree, which is 
equal to positive number (i + 2) (see e.g. P]) 



A' 



i + 2, i > 0. 



(69) 
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For instance, the component r]i of the basis solution r]\ of determining equations (^Tj) in 
accordance with (^6|), ( |67|) and (|69| ) has the homogeneity degree equal to four 



0\2 , 



4. 



(70) 



The indexing of the presented infinite (countable) vectors by one more integral number 
k > yields the desired representation of all linear independent solutions of determining 
equations (^) in the form of tensor of the second rank (matrix) ry^, in which the lower 
index A; > indicates the index of the basis infinitesimal operator in the general element 
of an infinite Lie algebra under consideration 



i,k=0 



(71) 



Under the conditions (|67| ) the integration of determining equations (^) for the given basis 
vector ?7| for a fixed value /c > is carried out with boundary conditions, that are imposed 
by requirements ( p7| ) in a single way. 

The representation of matrix 77^ for different lines are as follows {i is the column 
number, k is the line number) 



Vk 



{0, . . . , 0, 1, 0, -{k + 1)^0, -{k + 1)A\ ...}. 



(72) 



Here zeroes preceding unity describe matrix elements, which exist only for i < k, i.e. 
which are located below the principle diagonal i = k, that contains only units. The first 
nearest upper off-diagonal i = k + 1 also contains only zeroes. Expressions for elements 
from the second i = k + 2 and the third i = k + 3 upper off-diagonals are given in ( |72|) 
explicitly: they contain monomials, the homogeneity degree of which is equal to 2 and 3 
respectively, while the numerical coefficient {k + 1) is defined by the line number. 

In general, any one of the nonzero off-diagonals i = k + s with the number s > 2 is 
presented by polynomials with the homogeneity degree equal to s. This "line scheme" 



(72) is readily illustrated by a pictorial rendition of elements of the high left block of the 
discussed matrix (0 < i < 5, < /c < 3) 



( 1 





-^0 




^ (^0)2 


-A^ + 2A^A^ . 







1 





-2^0 


-2A^ 


-2^2^3(^10)2 . 










1 





-3^0 


-3^1 













1 





-4^0 





V 



(73) 



As a more illustrative example we present here the element ri\ of the matrix (|7^ with 
sufficiently high column number i = 10 and the homogeneity degree 9, that is located in 
the line with k = 1 (the second from above) 



^1 



10 



-2^7 + &A^A^ + &A^A^ + &A^A^ - 12^3(^0)2 
-2AA^A^A^ - A{A^f + 2{)A^{A^f . 



(74) 
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4.4 Illustrative example for matrix elements 

A much more comprehensive idea of definite expressions of matrix elements r/^ is given by 
the following list of elements (with the previous result included) of the first 11 columns 
(0 < i < 10) and 4 lines (0 < A; < 3) of matrix r/^, which define the k-th basic solution 
of determining equations (|6l| ) for vectors of the canonical infinitesimal operator (^6[), 
(p8|). The lower index "/c" is omitted for simplicity. 

0) A: = 0; r/° = 1, r]^ = 0, [rj'] = i, i > 2. 
rf = , 
rf" = -A^ , 

4 

r] 

rf 

rf = -A^ + 2A^A^ + 2A^A^ - 3^^ {A^f , (75) 

^8 ^ + 2^0^4 + 2^3^! + (^2)2 - 3^2(^0)2 

-3^0(^1)2 + (^0)4, 

7.9 = ^ 2A^A^ + 2^4^! + 2^3 ^2 - 3^3(^40)2 

- QA^A^A^ - {Ay + 4^1(^0)3 , 

7.10 = _^8 + 2A^A^ + 2A^A^ + A\2A^ - 3(^0)2] + A^\A^ - ^A^ A^\ 

+ ^2[_3(^l)2 _ 3^0^2 ^ 4(^0)3] ^ 6(^1)2(^0)2 _ (^0)5 _ 



.4 _ _^2 ^ (^0)2 ^ 
.5 _ ^ 2^0^! 



1) h=\\ rf = 0, ry^ = 1, rf = 0, [r?*] = i - 1, i > 3. 
rf = -2A^ , 
rf = -2A^ , 
rf = -2^2 ^ 3(^0)2 ^ 
.6 _ _2^3 ^ g^O^l ^ 

-2^4 + 6^0^2^ 3(^1)2 _ 4(^0)3^ (7Q) 



7?^ 
7 

v 

rf = -2A^ + QA^A^ + QA^A'' - 12^^(^0)2 

Q 



^10 



2A^ + GA^A"^ + 6^3^! + ^2 [3^2 _ i2(^0)2] 
12^0(^1)2 + 5(A0)\ 

-2^7 + Q^O^S + 6^4^! + 6^3 [^2 _ 2(^0)2] 

- 24A'^A^A'^ + ^1 [-4(^1)2 + 20(^0)3] . 



2) k = 2] rf = 0, r/i =0, ry2 = 1, 7^3 = g, [r/*] = i - 2, z > 4. 
= -3^0 ^ 

7^5 = -3^1 , 

7^6 = _3^2 ^ g(^0)2 ^ 

7^7 = -3^3 ^ 12^0^1 ^ (77) 

7^8 = -3^4 ^ 12^0^2 ^ g(^l)2 _ 10(^0)3 , 

rf = -3^5 ^ 12^0^3 ^ 12^2^1 _ 30^1(^0)2 ^ 

7^10 = -3^6 ^ 12^0^4 ^ 12^3^1 ^ g(^2)2 

- 30^0(^1)2 + 15(^0)4 - 30^2(^0)2 . 
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3) k = 3; 


rf = 0, ??i = 0, 


rf = 0, ry3 


= 1, r/^ = 0, [rf 


rf' = - 


-4^0, 






= - 


-4^1 , 






ri' = - 


-4^2 + 10(^0)2 , 






rf = - 


-4^3 _^ 20^0^1 , 






Q 

■q^ = - 


-4^4 + - 


f 10(^1)2 - 


20(^0)3 , 




-4y4S + 20^° ^3 


+ 20^2^1 _ 


-60^1(^0)2. 



5 Conclusion 

This paper presents a result of calculation of the infinite (countable) part of Lie point 
group admitted by the system of Benney equations - moment equations (0). In standard 
(non-canonical representation) the point Lie group of Benney equations (y) is described 
by the infinitesimal operator 

x=edt+ed^+Y.^'f^,. (79) 

where coordinates ^ and r/ obey the system of determining equations 

oo 
3=0 



(80) 



rjl+rjl+^+iA^-^r^o = 0, i, k > . 
Determining equations (pO|) result from (E^) in account of relationships between coordi- 



nates of infinitesimal operators (7£) and (pq) 

ae^ = V + eiK^' + iA'-'Al) - fA^ . (81) 



Infinitesimal operators (|34[), that were presented above, gives rise to the following coordi- 
nates 

(82) 

^ iA^-iK"^ + (i + 2)A'{K^ - K^) . 

The problem of finding coordinates of the operator (^) was first treated in |Q, where 
only these solutions, namely (^), ( |3^ and (^2[), were described. The main result of our 
paper is that point symmetries of Benney equations @) are exhausted by formulas ( |3^ 
and solutions of determining equations (|6ll), i.e. determining equations (^) do not have 
any other solutions. Solutions of determining equations ([GI] ) which are responsible for the 
infinite part of the point group probably have not been known so far. 
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As a next step it seems intriguing to generalize the result (^), i.e. to find the first 
order Lie-Backlund group admitted by Benney equations (^) with coordinates ae* of the 
canonical infinitesimal operator (|36[), that has the linear form 



8e* = ^' + E^'''4, ^>0. (83) 

j=0 



Though the unique existence of the linear form ( |83D as well as the complete solution of 
determining equations]^ for the tensor r/*'-' has not yet been obtained, all known facts are 
in agreement with this linear form. In particular, results of ^ mentioned above are 
consistent with the following expression for the tensor r/*'-' of the linear form 

CO s—j—2 

rf^^ = J2 kHXA+k,j+i + s J2 + k)A^^^-^HX^l^, ; i, j, s>Q. (84) 

k=Q k=0 

Here s is the number of the basis solution (similar to that used for r/* in (ff^)), is a 
polynomial of the homogeneity degree {s + 2) in moments A^. Compatibility conditions 
for determining equations for the tensor ij^'^ give rise to many relationships for H^, for 
example 

oo 

Y,jA^-'H%=sH^~\ s>0. (85) 

j=0 

An explicit form for the polynomial is presented below just to illustrate the aforesaid 

= A'^ + 7A^A^ + 7A^A^ + 7A^A^ + 21^3(^0)2 + 42^2^^° 

+ 7{A^f + 35A^{Ay. ^^^> 



Comparison between formulas (|7^) and ( p6D shows that they differ only in numerical values 
(and signs) of coefficients. The generating function for polynomials H'^ is given in m 
So constructing of a recursion operator, which transforms the linear form (^) for point 
group to the linear form (|8^ ) is of principal interest. 

To complete the conclusion we present formulas for the infinite dimensional algebra of 
a Lie point group admitted by the system of equations (P) 



Xi = + I 65x + I {xitt - vit) d.-^-itdg - \ (2{g + 2h)it + x\ut) dh , 

(87) 



2 2' - 2 y 4 

X2 = x{t)dx + Xtdv - xxttdh , ^3 = xdx + vdy + gdg + {g + 2h)dh 
X4 = fdf + g {dg -dh) , X5 = df + v [dg - dh) , Xq = n{t)dh , 

Xr = G{t,x,g + h){dg-dh) . 

Here in ( p7[ ) ^(t), x(0> /^(O G{t,x,g + h) are arbitrary functions of their arguments. 
Using the procedure described in section it is easily checked that prolongation of 

^For simplicity these equations are omitted here. 
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generators (jl^) on nonlocal variables g and h produce generators that directly follow from 
(|87|). Namely, the prolongation of Xi from (15) gives Xi provided ^ = 1, and X2 from (|l5|) 
gives X2 provided x = 1- Next, prolongation of X3 from (15) gives X2 provided x = *i 
and X4 from (15) gives Xi — X3/2 — X4 provided = t. At last, prolongation of X^ from 
(111) gives X3 + X4. 



It should be noted that on one hand the transition from nonlocal Benney equations 
(|lj) to their "potential" analog does not introduce additional independent variables that 
relate to restriction onto the boundary sets (compare with WE, 14|). This fact reflects the 



property of the fast decay of the distribution function / in Benney equations at the infinity, 
/(+CX)) = /(— 00) = 0. On the other hand the generators (^) give one more evidence in 
favor of constructing symmetry of nonlocal equations using various representations and 
different approaches. 



Acknowledgements 

This work is supported in part by a grant from the National Research Foundation (NRF) 
of South Africa and by the Department of Mathematics, Blekinge Institute of Technology. 
One of the authors (VFK) is also supported from the RFBR grant 99-01-00232, 00-15- 
96691. We thank our colleagues Sergey Meleshko and Rafail Gazizov for their valuable 
comments. 



References 

[1] Ibragimov, N.H., Elementary Lie Group Analysis and Ordinary Differential Equa- 
tions, John Wiley & Sons, Chichester, 1999. 

[2] Ibragimov, N.H., Introduction to Modern Group Analysis, Tay, Ufa (Russia), 2000. 

[3] Benney, D.J., 'Some properties of long nonlinear waves'. Studies in Applied Mathe- 
matics, Lll(l), 1973, 45-50. 

[4] Krasnoslobodtzev, A.V., 'Gas dynamic and kinetic analogies in the theory of vertically 
inhomogeneous shallow water'. Transactions of Institute of General Physics, USSR 
Academy of Sciences 18, 1989, 33-71 [in Russian]. 

[5] Kupershmidt, B.A. and Manin, Yu.L, 'Long- Wave Equation with Free Boundary. I. 
Conservation Laws and Solutions', Functional Analysis and its Applications 11(3), 
1977, 188-197. 

[6] Kupershmidt, B.A. and Manin, Yu.L, 'Long- Wave Equation with Free Boundary. II. 
Hamiltonian Structure and Higher Equations', Functional Analysis and its Applica- 
tions 12(1), 1978, 20-29. 

[7] Zakharov, V.E., 'Benney Equation and Quasiclassical Approximation in the method 
of the inverse problem', Functional analysis and its applications, 14(2), 1980, 89-98. 



19 



[8] Gibbons, J., 'Collisionless Boltzmann equations and integrable moment equations', 
Physica D3 3(3), 1981, 503-511. 

[9] Taranov, V.B., 'On the symmetry of high-frequency plasma motions (kinetic theory)'. 
Preprint KIYaI-74-21, Kiev, 1974. 

[10] Taranov, V.B., 'On the symmetry of one-dimensional high-frequency motions of a 
collisionless plasma', Journal of Technical Physics 46, 1976, 1271-1277 [in Russian]. 

[11] Taranov, V.B., 'Continuous symmetries of longitudinal motions of a collisionless 
plasma, invariant solutions and conservation laws'. Preprint of Institute of Theo- 
retical Physics, ITF- 78-161, Kiev, 1979 [in Russian]. 

[12] Bunimovich, A.I. and Krasnoslobodtsev, A.V., 'Group invariant solutions of kinetic 
equations'. Mechanics of fluids and gazes, No. 4, 1982, 135-140 [in Russian]. 

[13] Bunimovich, A.I. and Krasnoslobodtsev, A.V., 'On some invariant transformations 
of kinetic equations', Vestnik, Moscow State University, ser. 1., Matemat. mechan., 
No. 4, 1983, 69-72 [in Russian]. 

[14] Chctverikov, V.N. and Kudryavtsev, A.G., 'A Method for Computing Symmetries 
and Conservation Laws of Integro-Differential Equations', Acta Applicandae Mathe- 
maticae 41, 1995, 45-56. 

[15] Chetverikov, V.N. and Kudryavtsev, A.G., 'Modeling Integro-Differential Equations 
and a Method for Computing their Symmetries and Conservation Laws', American 
Mathematical Society Translations 167, 1995, 1-22. 

[16] Grigor'ev, Yu.N. and Meleshko, S.V., 'Investigation of invariant solutions of the non- 
linear Boltzmann kinetic equation and its models', Institute of theoretical and applied 
mckhanics, Siberian Division of Academy of Sciences, Preprint No. 18-86, Novosibirsk, 
1986. 

[17] Grigor'ev, Yu.N. and Meleshko, S.V., 'Group analysis of integro-diffcrential Boltz- 
mann equation', Doklady AN SSSR 297(2), 1987, 323-327 [in Russian]. 

[18] Meleshko, S.V., Doctor thesis. Novosibirsk, Institute of theoretical and applied 
mekhanics, Siberian Division of Academy of Sciences, 1991. 

[19] Meleshko, S.V., 'Application of group analysis in gas kinetics', Proc. Joint 
ISAMM/FRD Inter-disciplinary Workshop Symmetry analysis and mathematical 
modelling, 8-10 December 1998, N.H. Ibragimov and E.A. Lottering (eds.). Published 
by Foundation for Research Development (FRD, South Africa) and International In- 
stitute for Symmetry Analysis and Mathematical Modelling (ISAMM), 1998, 45-60. 

[20] Grigoryev, Yu.N., Meleshko, S.V., 'Group theoretical analysis of the kinetic Boltz- 
mann equation and its models'. Arch. Mech., 42(6), 1990, 693-701. 



20 



[21] Grigoriev, Yu.N., Meleshko, S.V., 'Bobylev-Krook-Wu Modes for Multicomponent 
Gas Mixtures', Phys. Rev. Lett., 81(1), 1998, 93-95. 

[22] Kovalev, V.F., Krivenko, S.V., and Pustovalov, V.V., 'Group symmetry of the kinetic 
equations of a collisionless plasma', JETF Letters 55(4), 1992, 256-259. 

[23] Kovalev, V.F., Krivenko, S.V., and Pustovalov, V.V., 'Group analysis of the Vlasov 
kinetic equation', Differential Equations 29(10), 1993, 1568-1578; 29(11), 1993, 
1712-1721. 

[24] Kovalev, V.F., Krivenko, S.V., and Pustovalov, V.V., 'Symmetry group of integro- 
differential equations', P.N.Lebedev Physical Institute Russian Academy of Sciences, 
Preprint No. 60, Moscow, 1998 [in Russian]. 

[25] CRC Handbook of Lie Group Analysis of Differential Equations. Volume 2. Appli- 
cations in engineering and physical sciences, N.H.Ibragimov (ed.), CRC Press, Boca 
Raton, FL, 1995, 396-431. 

[26] Kovalev, V.F., Krivenko, S.V., and Pustovalov, V.V., 'Symmetry group of Vlasov- 
Maxwell equations in plasma theory'. Journal of Nonlinear Mathematical Physics 
3(1-2), 1996, 175-180. 

[27] Gurevich, A.V. and Pitaevski, L.P., 'Nonlinear dynamics of a rarefied plasma and 
ionospheric aerodynamics'. Problems of Plasma Theory, Moscow, "Nauka", 1980, 
10, 3-87. 

[28] Bobylev, A.V. and Ibragimov, N.H., 'Relationships between the symmetry proper- 
ties of the equations of gas kinetics and hydrodynamics', J. Mathematical Modeling, 
1(3), 1989, 100-109 (English transl. in Mathematical Modeling and Computational 
Experiment, 1(3), 1993, 291-300). 



21 



